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1 Introduction 



In recent years, the Calogero- Sutherland type iV-body problems Jl|, in one 
dimension have attracted considerable attention not only because they are 
exactly solvable || but also due to their relationship with (l+l)-dimensional 
conformal field theory, random matrix theory || etc. In particular, the 
connections between exactly solvable models || and random matrix theory 
|| have been very fruitful. For example, by mapping these models to ran- 
dom matrices from an orthogonal, unitary or symplectic Gaussian ensemble, 
Sutherland was able to obtain all static correlation functions of the corre- 
sponding many body theory. The key point of this model is the pairwise long- 
range interaction among the N particles. One may add here that the family 
consisting of exactly solvable models, related to fully integrable systems, is 
quite small || and their importance lies in the fact that their small per- 
turbations describe wide range of physically interesting situations. Further, 
recent developments J7| relating equilibrium statistical mechanics to random 
matrix theory owing to non-integrability of dynamical systems has made the 
pursuit of unifying seemingly disparate ideas a very important theme. The 
results presented in this paper belong to the emerging intersection of several 
frontiers like quantum chaos, random matrix theory, many-body theory and 
equilibrium statistical mechanics ||. 

The universality in level correlations in linear (Gaussian) random ma- 
trix ensembles agrees very well with those in chaotic quantum systems || as 
also in many-body systems like nuclei ||. On the other hand, random matrix 
theory was connected to the world of exactly solvable models when the Brow- 
nian motion model was presented by Dyson ||10|| , and later on, by the works 
on level dynamics |Tl[. However, there are dynamical systems which are 



neither chaotic nor integrable - the so-called pseudo-integrable systems [12 



It is known that the spectral statistics of such systems are "non-universal 
with a universal trend" . In particular, for Aharonov-Bohm billiards, the 



level spacing distribution is linear for small spacing and it falls off exponen- 



tially for large spacing [14]. Similar features are numerically observed for the 



Anderson model in three dimensions at the metal-insulator transition point 



I5| . To understand these statistical features, and in the context of random 
banded matrices, a new random matrix model (which has been called as 
the short-range Dyson model in |T^|) was introduced |T7|, [TE| wherein the 



energy levels are treated as in the Coulomb gas model with the difference 
that only nearest neighbours interact. This new model explains features of 
intermediate statistics |16[ in some polygonal billiards. 

In view of all this it is worth enquiring if one can construct an iV-body 
problem which is exactly solvable and which is connected to the short-range 
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Dyson model (SRDM)? If possible, then using this correspondence one can 
hope to calculate the correlation functions of the corresponding many-body 
theory and see if the system exhibits long-range order and/or off-diagonal 
long-range order. 

The purpose of this paper is to present two such models in one dimension, 
one on a line and the other on a circle. We obtain the exact ground state 
and a part of the excitation spectrum on a line and the exact ground state 
on a circle in case the N particles are interacting via nearest and next-to- 
nearest neighbour interactions Further, in both the cases we show how 
the norm of the ground state wave function is related to the joint probability 
density function of the eigenvalues of short-range Dyson models. Using this 
mapping, we obtain one- and two-point functions of a related many-body 
theory in the thermodynamic limit and prove the absence of long-range order 
in the system. However, quite remarkably, we prove the existence of an off- 
diagonal long-range order in the symmetrized version of the corresponding 
many-body theory |20 . 



We also extend this work in several different directions. For example, 
we consider an iV-body problem with nearest and next-to-nearest neighbour 
interaction in an arbitrary number of dimensions D and show that the ground 
state and a part of the excitation spectrum can still be obtained analytically. 
We also obtain a part of the bound state spectrum in one dimension (both on 
a full line and on a circle) by replacing the root system An_i by BCntD^ 
etc. Besides, we also consider a model in two dimensions for which novel 
correlations are present in the ground as well as the excited states. 

The plan of the paper is the following. In Sec. II we consider an iV-body 
problem on a line characterized by the Hamiltonian (throughout this paper 
we shall use % = m = 1) 

i N o2 N-l I N-l -i 



2 dx^ (Xi (Xj_i Xj)(Xj Xj-i-i) 



with G > while g > —1/4 to prevent the collapse that a more attractive 
inversely quadratic potential would cause. We show that the ground state 
and at least a part of the excitation spectrum can be obtained if 

2 N 

g = P(P-l), G = p 2 ,V = ^-J2xl (2) 

A i=i 

Note that with the above restriction on G and g, (3 > 1/2. Further we also 
point out the connection between the norm of the ground state wave function 
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and the joint probability distribution function for eigenvalues in SRDM. In 
Sec. Ill we consider another iV-body problem, but this time on a circle 
characterized by the Hamiltonian 



H 



1 



N 







2§^ 2+ ^2gsin 2 [f(x J -x m )] 

2 N 



2 N 



- G j2 E COt 



< - \- 



cot 



(x N+ i = xi) , (3) 



(where again G > while g > —1/4) and obtain the exact ground state in 
case g and G are again as related by eq. (Q). Further, we also point out the 
connection between the norm of the ground state wave function and the joint 
probability distribution function for eigenvalues of short-range circular Dyson 
model (SRCDM). Using this connection, in Sees. IV and V we obtain several 
exact results for the corresponding many-body theory in the thermodynamic 
limit. In particular, in Sec. IV we calculate the two-particle correlation 
functions of a related many-body theory in the thermodynamic limit and 
prove the absence of long-range order in the system. In Sec. V we consider 
the symmetrized version of the model considered in Sec. Ill and show the 
existence of an off-diagonal long-range order in the bosonic system in the 
thermodynamic limit. In Sec. VI we consider the BCn generalization of 
the model (1) and obtain the exact ground state of the system. In Sec. VII 
we consider the BCjsr generalization of the model (3) and obtain the exact 
ground state of the system. In Sec. VIII we consider a generalization of the 
model (1) to higher dimensions and obtain the ground state and a part of 
the excitation spectrum. In Sec. IX we consider a variant of the model (HD in 
two dimensions and obtain the ground state as well a class of excited states 
all of which have a novel correlation built into them. Finally, in Sec. X we 
summarize the results obtained and point out several open questions. 



2 N-body problem in one dimension on a line 

Let us start from the Hamiltonian ([!]) and restrict our attention to the sector 
of configuration space corresponding to a definite ordering of the particles, 

say 



On using the ansatz 



1,2,..,JV-1. 



N-l 



$ = <P II ( X i ~ X i+1< 



(4) 



(5) 



8=1 
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in the corresponding Schrodinger equation Hip = Eip, it is easily shown that, 
provided g and G are related to (3 by eq. (@), satisfies the equation 

-2ga^-"g(i^)(^-a^J +(K -^ = ' (6) 

Following Calogero we start from <p as given by eq. @ and assume that 

= P,(x)$(r) . (7) 
where r 2 = x 2 . The function, $ satisfies the equation 

$"(r) + [JV + 2A; - 1 + 2(JV - l)/?]-$'(r) + 2[£ - V(r)]$(r) = , (8) 

provided Pfc(a;) is a homogeneous polynomial of degree k (k — 0, 1,2, ...) in 
the particle-coordinates and satisfies generalized Laplace equation 



N d 2 i /a d 



P fc (x) = 0. (9) 



We shall discuss few solutions of the Laplace equation (§) below. 

Let us now specialize to the case of the oscillator potential i.e. V(r) = 
\r 2 . In this case, eq. @ is the well known radial equation for the oscillator 
problem in more than one dimension and its solution is 

$(r) = exp(-Lur 2 /2)L*(ujr 2 ), n = 0, 1, 2, .... (10) 

where L^(x) is the associated Laguerre polynomial while the energy eigen- 
values are given by 



E, 



n 



N 

2n + k + — + (N - 



uj = Eq + (2n + k)u , (11) 



with a = — — 2n — 1. Few comments are in order at this stage. 
1. For large N, the energy E is proportional to iV so that 

lim ^ = (f3 + l)u, (12) 



N-+<x> N X 2 

i.e., the system has a good thermodynamic limit. In contrast, notice 
that the long-ranged Calogero model does not have a good thermody- 
namic limit since in that case for large N, E/N goes like N. 

2. The spectrum can be interpreted as due to noninteracting bosons (or 
fermions) plus (n, k)- independent (but iV-dependent) shift. 
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The ground state eigenvalue and eigenfunction of the model is thus given 
by (n = k = 0) 



E 



(N-l){1 + 



N 



CO 



(13) 



ipo = exp ^-7^ XX J II {xi-Xi+if . 



N 



N-l 



(14) 



A neat way of proving that we have indeed obtained the ground state can 
be given using the method of supersymmetric quantum mechanics [21|. To 
this end, we define the operators 



Qi 
Qi 

Qn 



d 



+ u)Xi + (3 



1 



2,3,..., N-l) 



d _ 1 
h UJXi - fj 

ax 1 x\ — X2 



d 



dx 



N 



+ LUX N + (3- 



(15) 



and their Hermitian conjugates Qf. It is easy to see that the Q's annihilate 
the ground state as given by eq. fll4|) . Further, the Hamiltonian ([[]) can be 
written in terms of these operators as 



N 



H — En 



z i=i 



(16) 



where Eq is as given by eq. (0). Now since the operator on the right hand 
side is nonnegative and annihilates the ground state wavefunction as given 
by eq. flllD , hence E as given by eq. (0) must be the ground state energy 
of the system. 

On rewriting ip in terms of a new variable 



one finds that the probability distribution for iV particles is given by 



(17) 



/ N \ N-l 

\ i=l ) i=l 



;is) 



where C is the normalization constant. We now observe that for (3 = 1,2,4, 
this ip 2 can be identified with the joint probability density function for the 
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eigenvalues of SRDM with Gaussian orthogonal, unitary or symplectic en- 
sembles respectively. We can therefore borrow the well-known results for 
these ensembles |]ISL |T7| and obtain exact results about a many-body theory 
defined in the limit, iV — > oo, uj —>■ 0, Nu = finite which defines the density 
of the system. For example, as N — > oo, the one-point function tends to a 
Gaussian for any (3 [|17j and is given by 



iV x 2 

i? 1 (x) = -_exp(-— ), (19) 

V27RX 2 2(T Z 

where a 2 = ^ +1 - > . Other results about the many-body theory will be dis- 
cussed in Sees. IV and V. 

Finally, let us discuss the polynomial solutions to the Laplace equation 
So far, we have been able to obtain solutions in the following cases: 

(i) k = 2,N > 2 (ii)k = 3,N > 3 (iii)k = A,N > 4 (iv)k = 5,N > 
5 (v)k = Q,N>Q. 

Besides we have also obtained solutions for k = 4, 5, 6 in case N = 3, 
and for k = 5, 6 in case N = 4. We find that for k > 3, the demand 
that there be no pole in Pk(x) alone does not require Pk(x) to be completely 
symmetrical polynomial. However, for k = 3, 4 and N = 3, 4 it turns out that 
solution to Laplace eq. (^) exists only if Pfc(x) is a completely symmetric 
polynomial. We suspect that this may be true in general. On assuming 
completely symmetric Pk{x) we find that in all the above cases we have a 
one-parameter family of solutions. In particular the various solutions are as 
follows (it is understood that the particle indices i, j, k, ... are always unequal 
unless mentioned otherwise). 

(i) k = 2,N > 2 

N N 

P k {x) = a^x 2 + b^2xiXj , (20) 

i=l i<j 

with (3 given by 

l3 - (N-l)( b -2a) - < 21 > 

(ii) k = 3,N > 3 

N N N 

P k (x) = a^xf + b x i x j + c x i x i x h , (22) 

i=l i,i=l i<j<k 

where c = 3(6 — a) and (3 is given by 

~ 3a+(iV - 1)6 (23) 
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(iii) k = 4, N > 4 



where 



N N N 

P k {x) = a ^ 4 + 6 x \ x j + c^x\x) 

i=l ij'=l i<j 

N N 
~\~d ^ ^ XjXfc -\- 6 ^ ] %i%j%k x l , 
ij<fe i<j<k<l 



e = 6(c - 2a) , d = 6 + 2c - 4a , 

(iV + 4)6 + 2(iV - 2)c - 4(iV - 2)a + 2(N - l)(2a + b - c)(3 



and (3 is given by 

6a + (N- l)c 







(JV — 1)(6 — 4a) ' 
(iv) fc = 5, iV > 5 



N N N N 

Pk(x) — a^x\ + b x t x i + c Y x i x "j + d Y x \ x i x k 

i=l i,j=l i,j=l i,j<k 

N N N 

+e y ' X^XjXk + / ^ ] X^XjX^Xi 4~ g ^ ' jXk x l x m i 

k,i<j i,j<k<l i<j<k<l<m 

where 

e = 5c — 5a — 36 , d = 6 + 2c — 5a , 
/ = 12c - 15a - 96 , g = 30(c - a - b) , 
(57V - 7)c - 3(iV - 4)6 - 5(iV - 2)a 
+ (iV-l)(5a + 36-2c)/3 = 0, 

and (3 is given by 

10a + (N- l)c 
^~ (7V-l)(6-5a) ' 

(iv) k = 6, N > 6 
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N N N N N 

Pk(x) = xf + b ^2 x i x j + c ^2 x t x ] + ^ x t x j x k + 

i=l i , j = 1 ij=l i,j<k i<j 

N N N N 

+/ Xj^XjXk -\- g XiXjXkXi -\- h x^XjXf.-\-p x^x^x^xi 

ij,k=l i,j<k<l i<j<k i<j,<k<l 

N N 
~\~Q ^ ^ 3^ XjXf t XiX m ~\- T ^ ^ XiXjX^XiX m X n , (30) 

i,j<k<l<m i<j<k<l<m<n 



where 



3e = 46 - 2c + 6a + / , d = 6 + 2c - 6a , a = 2/ + 2c - 46 - 6a , 

a = 2/ + 9a - 46 - c , p = 5/ + 18a - 86 - 6c , 

g = 6(2/ + 9a - 46 - 3c) , r = 30(/ + 6a - 26 - 2c) , 

(5N - 9)f - 2(4iV - 15)6 + 18(iV - 5)a - 6(iV - 5)c 

+(N - 1)(86 + 6c - 2f - 18a)/? = , 

146 - 2c + 6a + (N - l)f + 2(N - l)(3a + 26 - c)(3 = , (31) 
and (3 is given by 

It would be nice if one can find solutions for higher values of and further 
check if solutions exist (if at all) only if Pk{x) is a completely symmetric 
polynomial. While we are unable to prove it, we suspect that, subject to the 
solutions of the Laplace equation for higher k, we have obtained the complete 
spectrum for this problem. 

Finally it is worth enquiring if the bound state spectrum of the Hamilto- 
nian ([!]) can also be obtained in case the oscillator potential is replaced by 
any other potential. It turns out that as in the Calogero case |[22|| , in this 
case also the answer to the question is in affirmative. In particular, if instead 
the N particles are interacting via the iV-body potential as given by 

N l 

V(xi, x 2 , x N ) = -a^ I > (33) 

i=l \uli x i 



then also (most likely the entire) discrete spectrum can be obtained. This is 
because, after using the ansatz (^), eq. (Q) is essentially the radial Schrodinger 
equation for an attractive Coulomb potential and it is well known that the 
only two problems which are analytically solvable for all partial waves are 



10 



the Coulomb and the oscillator potentials. In particular the solution of (|8]) 
is then given by (note r 2 = Y,iLi x 1) 



$(r) = exp(-^2 | E \r)L» n {2p \ E \r) , (34) 
and the corresponding energy eigenvalues are 



E, 



a 2 



n,k 



n + k + f - 1 + (N - l)(3 



n 2 ' 



(35) 



when b = N+2k—3+2(N—l)/3. It may be noted that whereas in the oscillator 
case the spectrum is linear in (3, it is (— E)~ 1 / 2 which is linear in f3 in the 
case of the Coulomb-like potential. Secondly, as in any oscillator (Coulomb) 
problem, the energy depends on n and k only through the combination 2n + k 
(n + k). 

Is there any underlying reason why one is able to obtain the discrete 
spectrum for the iV-body problem with either the oscillator or the Coulomb- 
like potential ( |3~3"D ? Following p3| it is easily shown that in both the cases 
one can write down an underlying SU (1,1) algebra. Further, since the many- 
body potential W in ([[]) is a homogeneous function of the coordinates of 
degree -2, i.e. it satisfies 

N BW 

E^ = "2W, (36) 
i=i ux i 



hence, following the arguments of [23|, one can also establish a simple alge- 



braic relationship between the energy eigenstates of the iV-body problem ([I]) 
with the Coulomb-like potential ( |33D and the harmonic oscillator potential. 

It may be noted that the Hamiltonian (|I|) is not completely symmetric in 
the sense that whereas all other particles have two neighbours, particle 1 and 
iV have only one neighbour. Can one make it symmetric so that all particles 
will be treated on the same footing? One possible way is to add some extra 
terms in H. For example, consider 

H 1 = H + H', (37) 

where H is as given by eq. ([[]) while H' has the form 

1 -(38) 



H ' = 7 ~ — \2~ G 

(x N - x x y 



(x N - xi)(xi - x 2 ) + (sat_i - x N )(x N - Xi) 



Clearly, by adding these extra terms, the problem has become cyclic invariant 
for any iV while for N = 3 it is identical to the Calogero problem and hence 
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is in fact completely symmetric under the interchange of any two of the three 
particle coordinates. It may be noted that in the thermodynamic limit, these 
extra terms are irrelevant. 

We can again start from the ansatz (||) (but with N — l replaced by N) in 
the Schrodinger equation Hiip = Eip and using eq. @ we find that <fi again 
satisfies eq. (|6]) but with N — 1 in the second term being replaced by N. On 
further using the substitution as given by eq. (|7|) one finds that $ satisfies 
eq. (|8|) but with the coefficient of the 2(3 term being N instead of N — 1 
while Pk{x) is again a homogeneous polynomial of degree k (k=0,l,2,...) in 
the particle coordinates, which now satisfies instead of eq. (BJ) 



r N d 2 " 1 ( d d 



P k (x) = 0, (39) 



with xn+i = x\. 

How do the solutions to the Laplace eqs. |J) and ([39|) compare? For 
N = 3, eq. (|39"D is identical to that of Calogero and for this case Calogero 
has already obtained the solutions for any k. For iV > 3 and for k > 3, 
the demand that there be no pole in Pk{x) alone does not require Pk{x) to 
be completely symmetrical polynomial. However, for k = 3, 4 and N = 4 
it again turns out that solution to Laplace eq. (|39|) exists only if Pk{x) 
is a completely symmetric polynomial. We suspect that this may be true 
in general. On assuming completely symmetric Pk{x) we have been able to 
obtain a two-parameter family of solutions in case k=3,4,5,6 and N > k (note 
that for eq. @ we have obtained only one-parameter family of solutions). As 
an illustration, the solution for N > 4 and k = 4 is given by (it is understood 
that the particle indices i,j, k, ... are always unequal) 



N N N 

nr. . , 
1 3 



Pk(x) = x t + b ^ x\x,j + c^xfx 2 

i=l *,i=l i<3 

N N 

I d ^ x ^ x j x I c ^ XiXjXfcXi , (40) 

i,j<k i<j<k<l 

where 

e = 2(2a + 2d - 26 - c) , (41) 

6a + (N — l)c + p[8a - 26 + (2c - d)(N - 2)] = , (42) 

66+(iV-2)d+2/3[2(iV-l)6-2(iV-4)a + (A^-4)c-2(A^-l)d] = .(43) 

Solution to the new $ equation can be easily written down in case V(r) = 
or if it is given by eq. (|3"3"D. For example, it is easily checked that in the 



12 



former case the solution is again given by eq. (|TT| ) but the energy eigenvalues 
are now given by 

N 

E n:k = [2n + k + — + N(3}uj = E + (2n + k)u . (44) 

Similarly, in the later case, the solution is as given by eq. fl35|) except that 
in the term containing (5, N — 1 must be replaced by N. 

Apart from these two potentials, where we have obtained the entire bound 
state spectrum, there are several other potentials which are quasi-exactly 
solvable. For example, for the potential 



N / N \ 2 / N 

i=l \i=l / \i=l 



\ 3 
„2 \ 



(45) 



it is easily shown that the ground state energy and eigenf unctions are 



£ = _JL[]V + 2(iV-l)/3], (46) 



V> = exp 



n N r> N -, N-l 



, vE^) 2 + T7^E^ 2 Ilte-W, (47) 
4 j=1 4VC i=1 J i=1 

provided A, B, C are related by 

A= — - [AT + 2 + 2(iV-l)/5]v / C. (48) 

It is worth enquiring if the probability distribution for iV particles corre- 
sponding to (f47D can be mapped to some matrix model. In this context let 
us point out that the corresponding (long-ranged) Calogero problem was in 
fact mapped to the matrix model corresponding to branched polymers |24 . 



So far as we are aware of, answer to this question is not known in our case. 



3 N-body problem in one-dimension with pe- 
riodic boundary condition 

Soon after the seminal papers of Calogero JTJ and Sutherland where they 
considered an iV-body problem on full line, Sutherland [^5| also considered 
an iV-body problem with long-ranged interaction and with periodic bound- 
ary condition. He obtained the exact ground state energy and showed that 
the corresponding iV-particle probability density function is related to the 
random matrix with circular ensemble |2"5| . Using the known results for the 
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random matrix theory ||, he was able to obtain the static correlation func- 
tions of the corresponding many body theory. It is then natural to enquire 
if one can also obtain the exact ground state of an iV-body problem with 
nearest and next-to-nearest neighbour interaction with periodic boundary 
condition (PBC). Further, one would like to enquire if the corresponding N- 
particle probability density can be mapped to some known matrix model. 
The hope is that in this case one may be able to obtain the correlation func- 
tions of a related many-body theory in the thermodynamic limit. We now 
show that the answer to the question is in the affirmative. 

Let us start from the Hamiltonian (^]). We wish to find the ground state 
of the system subject to the periodic boundary condition (PBC) 



ip(x 1 , ...,Xi + L, ...,x N ) = ip(x 1 , ...,x u ...,x N ). 
For this, we start with a trial wave function of the form 



(49) 



N 

n su) 

i=\ 







-( 

L \Xi 



, (xtv+1 — X\) ■ 



(50) 



In this section, we restrict the coordinates Xi to the sector L > x\ > X2 > 
... > xn > 0, so that eq. (|50|) makes sense even for noninteger (3. The exten- 
sion to the full configuration space will be made in Sec. 5. On substituting 
eq. (|50|) in the Schrodinger equation for the Hamiltonian (§), we find that it 
is indeed a solution provided g and G are again related to (3 by eq. (0). The 
corresponding ground state energy turns out to be 



Np 2 n 2 

~1X 



En 



(51) 



The fact that this is indeed the ground state energy can be neatly proved 
by using the operators 



Qi 



d%4 L 



COt(Xi_i - Xi) - cot(^ - x i+l ) 



(52) 



and their Hermitian conjugates Q\ 
the ground state as given by eq. 
in terms of these operators as 



It is easy to see that the Q's annihilate 
. The Hamiltonian (§) can be rewritten 



H — En 



(53) 



where E is as given by eq. 
of the system. 



(|5*T|). Hence E must be the ground state energy 
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Thus unlike the Calogero-Sutherland type of models, our models (both 
of Sec. II and here) have good thermodynamic limit, i.e., the ground state 
energy per particle (= E /N) is finite as N — > oo. 

Having obtained the exact ground state, it is natural to enquire if the 
corresponding iV-particle probability density can be mapped to the joint 
probability distribution of some SRCDM so that we can obtain some exact 
results for the corresponding many-body theory. It turns out that indeed the 
square of the ground-state wave function is related to the joint probability 
distribution function for the SRCDM from where we conclude that the den- 
sity is a constant if < x < j-, and zero outside. Other exact results for the 
many-body theory will be discussed in the next two sections. 



4 Some exact results for the many-body prob- 
lem 

The square of the ground-state wavefunction of the many-body problem in- 
troduced in Sec. II (Sec. Ill) can be identified with the joint probability dis- 
tribution function of eigenvalues of the SRDM (SRCDM). Using SRCDM, 
Pandey |l7j and Bogomolny et al. |18| have shown that for any f3, the two- 
point correlation function has the form 

oo 

B^(s) = y £P^(n,8), (54) 

?1=1 

where s is the separation of two levels (or distance between two particles in 
the many-body theory considered here) and 

(R _|_ 1 \n(p+l) 
p(0)(rj A - {P + } nQ8+l)-l -Q8+]> 

F {U ' S) ~ r[n(/3 + l)] S 6 • (55j 

^From this expression it is not very easy to compute Ri{s) for arbitrary (3. 
However, it is easy to obtain the Laplace transform of R2{s) for any (3. In 
particular, if 



92(t) = / R2(s)e- ts ds, (56) 
J o 

then 

oo 

g 2 (t) = Y / 9(n,t), (57) 

n=l 

where g(n,t) is the Laplace transform of P(n, s), i.e., 

POO 

g(n,t) = / P(n, s)e~ ts ds . (58) 
Jo 
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On using P^(n, s) as given by eq. fl55| ) in eq. fl58|) it is easily shown that 



■ (59) 



(J3+l)n 



Hence 



1 



9f (t) = Y,9 m (n, t) = r!±m ; +1 _ , (60) 

n=l I 0+1 J 1 

from which one has to compute R^\s) by the Laplace inversion. 

For integer (3, it is possible to perform the Laplace inversion by making 
use of the fact that 

j j n— 1 ^2ikir/n 

17 I _ „2ifc7i7n ' (^1) 



x" — i n f— i x — e 2lkw / n 

fc=0 

yielding 

i?f ( S ) = f]fiV /3+1 ^- 1 ) (62) 

fc=0 

where 

ft = e 27rl/(/3+1) . (63) 

For (3=1, which corresponds to the orthogonal ensemble, the result is 
already known 0, |l| : R^\s) = 1 - e" 4s 



It is interesting to mention that R% (s) agrees very well with some of 
the pseudo-integrable billiards (e.g., the |-rhombus billiard). It is important 



here to note that for rhombus billiards | lf| , the Hamiltonian matrix has ele- 
ments which fall in their magnitude away from the principal diagonal. Thus, 
beyond a certain bandwidth, the elements are insignificant and the matrix is 
effectively banded. Immediately then, the results of banded matrices become 
applicable. Although there seems to be good agreement of the results from 



this random matrix theory as shown in [16|, []J|, in [lJJ it is also shown that 



there are other polygonal billiards for which R^ (s) is not an appropriate 
correlator. It is possible that for different bandwidths, and, by an inclusion 
of interactions beyond nearest neighbours in the short-range Dyson model, 
a family of random matrices result. This may, eventually, explain the entire 
family of systems exhibiting intermediate spectral statistics. 

Coming back to the two-point correlation function, depending on if (3 is 
an odd or an even integer, R2{s), as given by eq. (0), can be written in 
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a closed form which shows that B,2{s) is indeed real and further, it clearly 
exhibits oscillations for large s. In particular, it is easily shown that 



cos 



ran 



p+l 



+ (2p + 2)ssin 



m=l 

y P + T' 



(64) 



R 2 (P = 2p,s) = l + 2e-^ +1 >Y,e {2p+1)scos ^ ) 



m=l 



COS 



2mTT , . . 2m7r 

+ (2p+ l)ssm( 



2p+l w ' y 2p+ V 
For illustration, we give below explicit expressions for (5 = 2,3,4 



4 2) ( S ) 

4 3) (*) 

4 4) ( S ) 



7T 



1 — 2e 2 cos ■ . 

V 2 3j 3 



1-e 



-8s 



2e~ 4s sin(4s) 



1 + 2e 5s (- 1 + c -( 2 -/ 5 )) cos 

"47T 



2vr r . /2?r\ 

h 5s sin — 

5 V 5 J 



_l_ 2 e 5s ( _1 + cos ( 47r / 5 )) 



cos 



h os sin — 

5 V 5 



(65) 



(66) 



In Fig. 1, we have plotted R^\s) as a function of s for (3 = 1,2,3,4. 
These results show that, for integer j3, there is no long-range order in the 
corresponding many-body theory. 

Similarly, if (3 is half-integral, i.e., j3 = (2n + l)/2 then it is easily shown 
that 



Rp n+1)/2 \s) 



2n 



k=0 



2k -2zLtl s (i_o2fe) 



1 + erf 



(2n+l)s ^ fc 



, (67) 



where f2 is as given by eq. (|63|). 

For arbitrary /3, however, we are unable to perform the Laplace inversion 
and hence we do not have a closed expression for ^(s). However, one can 
numerically calculate it by using eqs. ([54] ) and (|55|). In Fig. 2, we have 
plotted R£\s) as a function of s for (3 = 1,4/3,3/2,5/3,2,7/3,5/2. From 
this figure it is clear that even for fractional /3, there is no long-range order. 
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5 Off-diagonal long-range order 



So far, nothing has been specified regarding the statistical character of the 
particles involved in the N-body problem of Sec. III. We now do that by first 
symmetrizing the Hamiltonian, that is by rewriting it as 

H = - X f ^ 

+ - ^P(2))-©(^P(JV-l) - XP{N))W{X P{1) , Xp (A ()^) 

PeS N 

where 6 is the step function and W(xi, xjy) is the N-body potential of eq. 
(3). Next, relying on the solution given in eq. (|50|), we introduce the (not 
normalized) wave function: 

iPn(xi, —,xn) = £p4>n(xp(i), —, Xp(N)) , (69) 

where P is the permutation in Sn such that 1 > xp(i) > £p(2) > ••• > £p(jv) > 
0,Ep = l(ep = sign(P)) in the N-boson (N-fermion) case and 

N 

<f> N {xi, x N ) = n | sinn(x n - x n+1 ) f ; (x N+1 = Xi) , (70) 

n=l 

(we have set the scale factor L equal to 1). Primitively, the function (|69|) is 
defined on the hypercube [0,1]^. The following properties of ipx are easily 
verified, provided that (3 > 2: 

1. In the bosonic case, ipN can be continued to a multi-periodic function 
in the whole space 1Z N (or equivalently on the torus T N ): 

^ N {xx, ...,Xi+l, ...,x N ) = ij} N {xx, -;Xi, -;X N ) ; (i = 1, N) ,(71) 

which belongs to C 2 (i.e. is twice continuously differentiable). Owing 
to this property and the results of Sec. 3, ipN then obeys the Schrodinger 
equation (with Hamiltonian (3) and energy as given by eq. (|5l|) ) not 
only in the sector X\ > x 2 > ... > xjq but every where. Thus, tjj^ 
describes the ground state wave function of the N-boson system. More- 
over, it is translation invariant (on 1Z ): 

i>N(%i +a,x 2 +a, ...,x N + a) = i/j n (x 1 ,x 2 , ■■■,x n ) ; V a e TZ. (72) 
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2. In the fermionic case, the continuation by periodicity is possible only 
for odd N, in which case eq. (|7l"[) still holds with ip^ e C 2 . For even N 
on the contrary, enforcing the periodicity ([71]) leads to a discontinuous 
function ip^, so that the Schrodinger equation is no longer satisfied on 
the configuration space T N . 

Therefore, in the following we shall implicitly restrict ourselves to odd 
values of N when treating fermions. The translation invariance (|72|) 
then remains valid. 

We are interested in the one-particle reduced density matrix, given by 
Pn(x-x') = — dxi... dx N -iip N (xi 7 ...,x N -i,x)ip N (xi,...,x N -i,x) 



C N Jo 

where CV stands for the squared norm of the wave function: 



C 



[ dx x ... [ dx N \ip N {x 1 ,...,x N ) | 2 . (74) 
Jo Jo 

That the R.H.S. of eq. ( [75] ) defines a (periodic) function of (x — x') is an 
easy consequence of eqs. (|7T|) and (|72|). The normalization of pjy is such that 
Pjv(O) = N, the particle density. Further, the function Pn{0 is manifestly of 
positive type on the U(l) group, which implies that its Fourier coefficients, 

fdZe-^pN®; (n = 0,±l,±2,...), (75) 
Jo 



(n) 

Pn 



are non-negative (Bochner's theorem). In fact, this directly appears if one 
writes their explicit expression 

iV r 1 , r 1 



(n) ^ f f 

Pn = ~FT I dxi... I dx N -iifj N (xi,...,x N -i,0)X 
Gat jo Jo 

X [ dxe 2innx ij N (x 1 ,...,x N _ 1 ,x), (76) 



in the form (obtained by using the periodicity property): 

(n) N f 1 f 1 f 1 

Pn = ~ft i dx\... / dx N -i \ / dxe 2mnx ip N (xi, xjv-i, x) \ 2 . (77) 
Cn Jo Jo Jo 

In the bosonic case, since the function is not only of positive type but 
also positive (like iPn), eq. ( [75] ) shows us that 

p!?>p!?; (n=±l,±2,...). (78) 

In the fermionic case, eq. (|78|) is not necessarily true (because ip^ changes 
sign on T ) and it is not an easy matter to determine the largest Fourier 
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coefficient. Notice that the coefficients p^\ which physically represent the 
expectation values of the number of particles having momentum k n = 2im in 
the ground state, are nothing but the eigenvalues of the one-particle reduced 
density matrix (diagonal in the k n representation). According to the Onsager- 
Penrose criterion |2TS|] , no condensation can occur in the system (at least for 
Bose particles) if the largest of these eigenvalues is not an extensive quantity 
in the thermodynamic limit, that is, if 

(o) 

li m PR- = o . (79) 

N^oo N 

For Fermi particles, this criterion is not sufficient, and one has to look also at 



the largest eigenvalue of the two-particle reduced density matrix [28|. Since 
we are presently unable to determine the largest eigenvalue of p^ itself in the 
fermionic case, we shall not discuss extensively the latter here. Nevertheless, 
we shall look for the large N behaviour of p$ for bosons and fermions at a 
time, as this does not require much extra work and can give some indications 
in the fermionic case too. Let us write: 



SO) 



N C N ' 
where Cat is given by eq. (|74[) and 

:/ dx x ... I dx N -ii) N (x 1 ,...,x N ^ 1 ,0)[ dxif) N (xi, ...,x N - U x) ,(81) 
Jo Jo Jo 



A 



(the expression (|76|) of p N ' is more convenient than (|77f ) for our purpose). 
Because of the special form (|69"1)-(|T0"|) of the wave function, the computation 
of the squared norm is already not a trivial task, in sharp contrast to 
the case of N free, impenetrable particles. As a consequence, the (mainly 



algebraic) method introduced long ago by Lenard [|27| to deal with the latter 
case does not apply here, and we have to resort to another device. For 
conciseness, we introduce the notation: 

S(x n - x n _i)= | sin7r(x n - x n+1 ) \ p , (82) 

and define: 

S 2 (A) = [ A dxS(x)S(A -x); (0 < A < 1) . (83) 
Jo 

Our starting point will be the following representations of Cat and An'- 



C N = (N- 1)!— / dxe~ ix F(x) N , (84) 

Z7T J-oo 
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A N = (N — 1)!— /°° dxe~ lx F(x 

ZTT J-oo 



(85) 



where 



and 



F{x) 
G{x) 
H(x) 

Vn = 



dAe iAx S(A) 



dAe iAx S 2 (A) 2 



dAe iAx S{A)S 2 (Af 



(N — 2) for bosons 
— 1 for fermions . 



(86) 
(87) 



The representations (p4|)-(p7|) follow from the convolution structure of the 
expressions (|73 ) and (|8T| ) of CV and A^, when written in terms of appropriate 
variables. Their proof is given in the Appendix. Our aim is to extract from 
them the large N behaviour of Cat and A^. Their form is especially suited 
for that purpose, because the integrands in eqs. (|84"D and (]85[) are entire 
functions, as polynomial combinations of Fourier transforms of functions with 
compact support (eq. (|8~6"D). Indeed, we are then allowed to, first, shift the 
integration path and then apply the residue theorem to meromorphic pieces 
of the integrands. However, it turns out that the calculations needed for 
arbitrary (integer) values of (3 are quite cumbersome. So, in order to keep 
the argument clear enough, we shall content ourselves to present below these 
calculations in the simplest case, namely (3 = 1 (recall that, strictly speaking, 
this value is not allowed), being understood that similar results are obtained 
for all integers (3 > 2. 

For (3 = 1, S(A) = shi7rA, and eq. (BO) gives, after reductions: 



F(x) 



G{x) 



Hix) 



2tt 2 1- 




i x(x 2 - 


- 4vr 2 ) ' 


4vr 4 5x 2 


-4vr 2 


i x 3 ix 2 


-4vr 2 ) 3 


4tt 3 


1 



+ e ix R { ~ l \ 



x 



% xix 2 — 47T 



2^2 



c R { ~ 2 \x) 



where R^i 
when x — > 

portance. This produces, for the functions to be integrated in eqs. ([84]) and 



x) is a generic notation for rational functions behaving like x n 
oo, and the precise form of which will be eventually of no im- 



Fix 



,2tt 2 



N 



^X 2 - 47T 2 )] 



N 



+ Y,e inx R ( n 3N) > 



X , 



n=l 



(89) 



21 



F(xf- 3 [F(x)G(x) + VnH(x) 2 ] 

Let us stress again that these functions, when analytically continued, are 
holomorphic in the whole complex plane (the poles appearing in the first 
term are exactly canceled by the remaining ones). 
We consider first Cn, now given by 

Since the function within the curly bracket is an entire one, we can shift the 
integration path to / = {z = x + ia \ x e TZ}. Let us choose a > 0. Then, 
by Cauchy theorem 

N 

dze~ iz J2 e mz R { - 3N) {z)=0. (92) 

n=l 

Indeed, the integrand is holomorphic above I and is bounded there by const. 
| z \~ 3N , which allows us to close the integration path at infinity in the upper 
complex plane. We end up with 



1 



Cn = (N-1)\[ — ) -\dz „, 2 _, 2W - (93) 



2ti z \ l_ 

2^ Ji"'z N (z 2 -4tt\ 

Similarly, we are allowed to close the integration path at infinity in eq. ( p3|) , 
but this time in the lower complex plane. The integrand has now poles at 
z = 0, ±27r, and applying the residue theorem leads to explicit expressions 
for Cat. Unfortunately, these expressions turn out to appear as (finite) sums 
with alternating signs, the terms of which become very close to each other for 
large N. They are therefore useless for determining the asymptotic behaviour 
of Cat, and we have to proceed differently. Let us write 

f a 7 - 1 rfJV-1 I A 7 f d~ e ~" 

•> I z N (z 2 —Att 2 ) n — (N-l)\da N -^ la=47r 2 J I az z N (z 2 -a) 

= (n%S-i U=4- 2 [R+(<x) + R-(<*) + R (a)} , (94) 

where R±(a) and Ro(a) are the residues of the last integrand at z = ±^/a 
and z = respectively. They are readily computed, assuming first that 
N = 2M + 1 is odd: 

t-, i \ n / n COSa/o ^ (-l) r r _ M _ ] 

R + (a) + R-(a) = -^ = £l_J_ a ' \ 

M (_-l\r 

Ma) = -E^"-"- 1 - (95) 
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Hence 



R+(a) + R-(a) + R (a 
Using eqs. (gD, (f~ 



1 



,M+1 



E 



c 



2tt 



" Q (2M + 2s + 2)! 
and (|96|) we then obtain 



ft' 



(96) 



N 



(N + n-l)\ 
f n\(3N + 2n-l)\ [ 



oo 

|q=47T 2 ^ t 



rft 



-47T 



s=0 (2M + 2. + 2)! 

2« ( g? ) 



The result is exactly the same for even N. It suffices now to observe that 
the last series alternates in sign and is decreasing to deduce 



(9* 



Our procedure for evaluating is quite similar, and we give below only 
the main steps. From eqs. fl85|) and fl90|) we get 



.4 



N 



{N-m^ri-Stdze-i 



1 (2tt 2 \N i d N 



N 



(2tt± \N j_ 

V i > 27T 



=47r 2 J 7 dze" 



: (5+2 TO )^ 2 -4tt 2 
Z n + 1 (z 2 -4tt 2 ) n + 1 

5+2r?jy 4tt 2 



(99) 



and, after computing the residues at z — ±y/a and z = 0, we get 



.4 



N 



(-2tt 2 ) n d N I 
N da N l Q = 



47T 2 



OO 

s=0 



4tt 2 



(V+2s)! (Af+2s+2)! 



-a)' 



(2^ 
AT 



n=0 



(Af+ra)! 



5+2t?iv 



4tH 



-47T 



2\n 



(3W+2n)! (3Ar+2n+2)! 

Again, the last series alternates in sign and decreases, which entails 

' N 

Finally, using eqs. (fO), ( 



Aat = (5 + 2r ?jV )(27r 2 ) jV v ';_ T ;, / ' [l + O 



(3N)\ 

, ( F01]) and fl87|) we obtain 



(100) 



:ior 



(o) 
Pn 

N 



5 + 2r] N 
3N 



|[l + 0(^)] for bosons 
4:[l + 0(4:)] for fermions 



(102) 



The same procedure applies for all integer values of /3, although the algebra 
becomes quite involved. The general result for bosons (and for any integer 
(3) is: 

pP _ (/3!) 4 [(3/3 + l)!] 2 



lim r 



[(2/?)!p[(2/3 + l)!]3 



(103) 
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Our method does not adapt straight forwardly to the case of non-integer 
values of (3, but there is clearly no reason to expect a different outcome for 
such intermediate values. Therefore, the Onsager-Penrose criterion ( 73 ) is 
not met for bosons, and we reach the conclusion that Bose-Einstein conden- 
sation is possible in the bosonic version of the iV-body model discussed in 
Sect. 3. 

In the fermionic version, the result (|102|) is not conclusive, as explained 
after eq. (^). It only points (not too surprisingly) to the absence of quantum 
phase in the system. 



6 The Bn model in one dimension 



Subsequent to the seminal work of Calogero and Sutherland for the A^-i 
system, the entire bound state spectrum of the Calogero model was obtained 
for BCn-i Dn root systems [§]. It is then natural to enquire if in our case, can 
one at least obtain the exact ground state and radial excitation spectrum in 
the BCn ot Dn case? We now show that the answer to this question is in 
the affirmative. 

Consider the BCn Hamiltonian, 



N 



H 



d 2 



N 



JV-1 r 



\i=l 



i=l 



1 



+ 



1 



iV-1 r 



i=2 
N 



1 



1 



(Xj — Xi + i) 2 (Xj + Xj + i) 2 

1 1 N 



+ 



i=l X i 



(104) 



of which B n ,Cn and D N are the special cases. We again restrict our atten- 
tion to the sector of configuration space corresponding to a definite ordering 
of the particles as given by eq. (|j). 
We start with the ansatz 



A' 



v>=^)0(r) n< 



x 



2W/2 



JV-1 



II K 2 - 



X 



i+lv 



(105) 



\i=l 



i=l 



where r 2 = Ya=iX 2 . On substituting it in the Schrodinger equation for the 
-BAT-Hamiltonian (1104!) we find that satisfies 



$ (r) + [iV+4£;-l+2iV7+4(iV-l) ( 5]-$ / (r)+2 [E - V(r)] $(r) = ,(106) 
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provided g and G are again related to (3 by eq. 
by 



while g\ is related to 7 



9i 



|(7-D. 



Here, P2k(x) is a homogeneous polynomial of degree 2k [k = 0, 1, 2, 
particle-coordinates and satisfies the generalized Laplace equation 



N q2 

Y — 

^ dx 

i=i Udj 



N 1 d 



7V-1 



-4/3 E 



x 



i+ly 



_d_ 

dxi 







dx 



(107) 
in the 



P2k(x) 



Let us now specialize to the case of the oscillator potential, i.e., V(r) = 
" 2 . In this case, ( |106| ) is the well known radial equation for the oscillator 



problem in more than one dimension and its solution is 
$(r) = exp(-ujr 2 /2)L*(ujr 2 ), n = 0, 1, 2, .... 



(109) 



where L^(x) is the associated Laguerre polynomial while the energy eigen- 
values are given by 



En 



N 

2n + 2k + — + N-f + 2(N- l)(3 



U) 



(110) 



with a = — — 2n — 1. The exact ground state is obtained from here when 
n = k = 0. The fact that n = k = gives the exact ground state energy of the 
system can be easily shown a la A^-i case by the method of supersymmetric 
quantum mechanics. It may be noted that for large N, the energy E is 
proportional to N so that like the A^-i case, the model also has a good 
thermodynamic limit. In contrast, notice that the long-ranged E>n Calogero 
model does not have a good thermodynamic limit. 

Are there homogeneous polynomial solutions of eq. ( |108| ) of degree 2k 
(k > 1)? While we are unable to answer this question for any k, at least 
for small values of k (k > 0) there does not seem to be any solution to eq. 
(|108|) . For example, we have failed to find any polynomial solution of degree 
2,4 and 6. Thus it appears that unlike the A^-i case, in the BCn case one 
is only able to obtain the ground state and radial excitations over it. 

Proceeding in the same way, the energy eigenvalues and eigenf unctions in 
the case of the Coulomb-like potential (|33) are 



E 



a 



n + 2k + ^ + N-f + 2(N- l)(3 



0.(108) 



$ = e-V^L b n {2^2 I E |r) 



(112) 



25 



where b = N — 2 + Ak + 2N / ~y + 4(N — l)/3. Again, so far we have been able 
to obtain solutions only in case k=0. 

As in Sec. II, in the BC^ Hamiltonian (|104[ ), all the particles are not being 
treated on the same footing. Again, one possibility is to add extra terms. 
Consider for example, 



H 1 = H + H', 
where H is as given by eq. ( |104| ) while H' has the form 

H' 



(113) 



1 



9 



G 



(x N - x{) 2 

, 1 

Xjq — X\ 



+ 



1 



(x N + Xi) 2 

1 )( 



1 



+ 



Xj\f + X\ 

1 



X\ — x 2 



)(- 



X\ 

+ 



2+ 



Xn-1 — Xn Xn-I+Xn Xn — Xl Xn + Xi 



(114) 



One can now run through the arguments as given above and show that the 
eigenstates for both the oscillator and Coulomb-like potentials have the same 
form as given above except that in the term multiplying (3, N—l gets replaced 
by iV at all places including in the Laplace eq. ( |108| ). However, now we find 
that there are indeed solutions to the Laplace eq. ( |108| ) (with N — l replaced 
by N). In particular, the solution for any iV(> 4) and k = 4 is given by 



N 



N 



Pk- 



=1 



where 



-2 



3 + 8/3 + 27 



N - 1 + 2(N - 1) 7 + 4(iV - 2)(3 



(115) 



;ii6) 



As in the A^-i case, we again find that even though the Laplace eq. ( |108| ) 
is only invariant under cyclic permutations, the solution is in fact invariant 
under the permutation of any two coordinates. It will be interesting to try 
to find solutions for higher values of k and study the full degeneracy of the 
spectrum. 

Besides these two, one can obtain a part of the spectra including the 
ground state for several other potentials but we shall not discuss them here. 
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7 BCn model in one dimension with periodic 
boundary condition 

Following the work of Sutherland [[25[ on the An-i root system, the exact 
ground state as well as the excitation spectrum was also obtained in the case 
of the BCn, Dm root systems ||. It is then worth enquiring if, in our case, 
one can also obtain the ground state and the excitation spectrum. As a first 
step in that direction, we shall obtain the exact ground state of the BCn 
model with periodic boundary condition. 

The Hamiltonian for the BCn case is given by 



^2 9 T 2 

1=1 



H 



2 N 

2^ ax? 1 ^ 
2 i 



i 



i 



sin 2 j(xi 



X i+1, 



sin 2 f(xi + x i+ i) 



7T _ 

L z , sin 4-x 



L ^ V 1 



7T 2 N 
L i=l 



7T 



Cot (^i — 1 



7T 



- COt — (Xi-i + Xi 



7T 7T 
COt —{Xi - X i+ i) + COt —{Xi + 



(117) 



We again restrict our attention to the sector of the configuration space cor- 
responding to a definite ordering of the particles as given by eq. @. For this 
case, we start with a trial wave function of the form 



N 



N 



N 



N 



^o = U shv' Oi II(sm 2 26i)^ 2 l[[sm 2 (e i -9 i+1 )f 2 n[sm 2 (^+ 1 )] /32 ,(118) 



i=l 



i=l 



i=l 



i=l 



{6i = irxi/L) and substitute it in the Schrodinger equation for the Hamilto- 
nian (|117f) . We find that it is indeed a solution provided g and G are again 
related to (3 by eq. (0) while 01,(72 are related to 7,71 by 



9\ 



7 



[7 + 271-I], & = 271(71-1) 



;ii9) 



The corresponding ground state energy turns out to be 
A% 2 



Eq 



2L 2 



( 7 + 7i + 2/?) 2 . 



(120) 



The fact that this is indeed the ground state energy can be easily proved as 
in Sees. II and III. 
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8 TV-body problem in D- dimensions 

Having obtained some results for the iV-body problem ([]]) in one dimension, 
we study generalization to higher dimensions. Let us consider the following 
model in D-dimensions : 



N 



N-l 



H 



" i=l 

N-l 



i=2 



i= i K n-r i+ iy 
(fj_i - ri).(ri - f i+1 ) 
(fj_i - fi) 2 (ri - r i+1 ) 2 



( N 

We 



:i2V 



\i=l 



On using the ansatz, 



N-l 

n 

i=l 



n - r 1+ i 



A? 



r , r 



E r 



2 ) 



122) 



in the Schrodinger equation for the Hamiltonian (121), it can be shown that 
4>{r) satisfies 



cf) (r) + [DN -1 + 2(N- l)0\-(j>'(r) + 2(E - V(r)U(r) = , (123) 

r 



provided g and G are related to (3 by 
g = f3 2 + (D-2)(3, G = f3 2 . 



(124) 

Equation (|123|) is easily solved in the case of the oscillator potential (i.e., 

V(r) 



r 2 ) yielding the energy eigenstates as 



UJ 



exp ( -^-r 2 ) L b n (ujr 2 



Er, 



2n + (N — 1)(3 + 



DN 



UJ . 



(125) 
(126) 



Here 6= £-2n-l. It may be noted that as in all other higher dimensional 
many-body problems, one has only obtained a part of the energy eigenvalue 
spectrum which however includes the ground state. In particular, the ground 
state energy eigenvalue and eigenfunction is given by 



E 



(N-1)P + 



DN 



(127) 



4>o = exp 



to 



N 



E-n n 



i=l 



N-l 



i=l 



n - r l+ i 
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As expected, for D = 1 these results go over to those obtained in Sec. II. 
The fact that this is indeed the ground state energy can be easily proved by 
using again a supersymmetric formulation |21[j . 

At this point it is worth asking if the probability distribution for N parti- 
cles (at least for some D(> 1)) can be mapped to some known random matrix 
ensemble ? In this context we recall that in the case of the Calogero-type 
model, it has been shown that in two space dimensions | ipo 1 2 can be mapped 



to complex random matrix |JT|. Using this identification one was able to cal- 
culate all the correlation functions of the corresponding many-body theory 
and show the absence of long-range order but the presence of an off-diagonal 
long-range order in that theory. Unfortunately, so far as we are aware of, 
answer to this question is unknown in this particular case. We hope that at 
least in the case of two space dimensions, where | ipo | 2 for our model is given 
by 




| ipo 1 2 can be mapped to some variant of the short-range Dyson model. 

Finally, we observe that the ground state and a class of excited states can 
also be obtained in D-dimensions in case the oscillator potential is replaced by 
the iV-body Coulomb-like potential V(r) = —a/ yJ2 r i, because the resulting 
equation ( |123| ) is essentially the radial equation for the Coulomb potential. 
In particular, the energy eigenvalues and eigenfunctions are given by 

En = " T2 , (130) 

2 n +^ + (N-l)P 

AT-1 



^ n = e W {-J2\E\r)L h n {2p\E\r)[ J[ \r i -r i+1 \n, (131) 



i=l 



where h' = DN — 2 + 2(iV — 1)(3. It may again be noted that whereas the 
ground state energy is linear in /3 in the oscillator case, it is not so in the 
case of the Coulomb-like iV-body potential. 

9 Short-range model in two dimensions with 
novel correlations 



Few years back, Murthy et al. |32 considered a model in two dimensions 



with two-body and three-body long-ranged interactions and obtained the 



29 



exact ground state and a class of excited states. The interesting feature of 
this model was that all these states had a built-in novel correlation of the 
form I Xaa \ 9 where 



(132) 



It is then natural to enquire if one can construct a model in two dimensions 
and obtain ground and few excited states of the system all of which would 
have a built-in short-range correlation of the form 



X jtj+ i — xjy j+ i - VjX j+ i 



(133) 



We now show that this is indeed possible. Let us consider the following 
Hamiltonian 



H 



1 



N 



2 N 



N-1 fS ,=S 
A i,i+1 



JV-1 



Tj-i ■ r i+ i 

,_2 ^i-l,i-^i,i+l 



2 ^ 

* i=l " i=l i=l 

where X ii+ i is as given by eq. ( 133 ). We start with the ansatz 



(134) 



N-1 



II X i,i+1 



(135) 



On substituting the ansatz in the Schrodinger equation Hip = Eip, one finds 
that satisfies the equation 



i iV jv xv— i -r 



N 



N-1 



i=l 



d 



i=l 

d 



dx 



Xi 



i+1 



1=1 

—)]<f>=(E-[N + 2(N-l)0\u)(f>, (136) 



x„ 



i,i+l 



Xi+l 



_d_ 

dyi 



Vi+i 



_d_ 

dxi 



provided g and G are related by (2). It is interesting to note that even 
though we are considering the novel correlation model in two- dimensions, 
the relationship between g and G is as in the case of our one-dimensional 
model. We do not know if this has any deep significance. 

We conclude from here that ip, as given by eq. (|135|) , with <p being a 
constant is the ground state of the system with the corresponding ground 
state energy being 



E = [N + 2(iV - 1) 



OJ. 



(137) 



Let us remark that, like the relationship between coupling constants, the 
ground state energy too has essentially the same form as that of the one- 
dimensional short-range A^-x model as given by eq. (|13D. That one has 
indeed obtained the ground state can be proved as before. 
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As in other many-body problems in two and higher dimensions, we are 
unable to find the complete excited-state spectrum. However, a class of 
excited states can be obtained from ( |136| ). To that end we introduce the 
complex coordinates 

. * . „ d 1/9 . 9 \ d 1/9 . 

z = z = x-iy, 9 = — = - - ^ ,9 = — = - — + z 

cte 2 ycte ay y c>z* 2 \ox 

In terms of these coordinates, the differential eq. (|136|) takes the form 



n N-i z i+1 di - Zid i+ i + z*d* +1 - z* +1 d* 

2j2d t d* + 2/3 X: 



i=l i=l 
A" 



Z i Z i+l 



0. (139) 



Now it is readily proved shown that the Hamiltonian H commutes with the 
total angular momentum operator L = Y^iLiizA — z*d*), so that one can 
classify solutions according to their angular momentum: Lip = 1(f). 

On defining t — u J2i z i z \ an d let <p = <p{t) it is easily shown that cj>{t) 
satisfies 



t<p"(t) + 



UJ 



m+(^-^)m=*> (140) 



where Eq is as given by eq. ( |137|) . Hence the allowed solutions with I = 



arc 



E = E + 2nuj, <j)(t) = L n " (t) . (141) 

Solutions with angular momentum I > or I < can similarly be obtained 
by introducing t z = ujJ2i z i or t z * = ^J2i( z i) 2 - F° r example, let (p = 4>(t z ). 
Then eq. (|139j) reduces to 



2ujt z -^ = (E - E )<f) . (142) 
dt z 

This is the well known Euler equation whose solutions are just monomials 
in t z . The solution is given by <f>(t x ) = t™(m > 0), and hence the angular 
momentum I = 2m while the energy eigenvalues are E = Eq + 2muj = 
Eq + Ioj. Further, we can combine these solutions with the I = solutions 
obtained above and obtain a tower of excited states. For example, let us 
define <p(zi,z*) = 0!(t)0 2 (t 2 ), where <p\ is a solution with I = 0, while $2 is 
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the solution with I > 0. On using (p2(t z ) = t™ it is easily shown that (f>i again 
satisfies a confluent hypergeometric equation, 



h 2m 

UJ 



t 



E-E 
2uj 



+ rnj 0i(t) = 0. 



(143) 



Hence the energy eigenvalues are given by E — Eq = (2n r + 2m)w. One 
may repeat the procedure to obtain exact solutions for a tower of states with 
I < 0. 



10 Summary 

In this paper we have discussed an iV-body problem in one dimension and pre- 
sented its exact ground state on a circle and most likely the entire spectrum 
on a real line. There are several similarities as well as differences between the 
model discussed here and Calogero-Sutherland (CSM) type of models and it 
might be worthwhile to compare the salient features of the two. 

1. Whereas in CSM the interaction is between all neighbours, in our case 
the interaction is only between nearest and next-to-nearest neighbours. 
Note however that in both the cases it is an inverse square interaction. 

2. Whereas in CSM (in one dimension) there is only two-body interaction, 
both two- and three-body interactions are required in our model for 
partial (or possibly exact) solvability on a real line. 

3. Whereas the complete bound state spectrum is obtained in the Suther- 
land model (periodic potential) or if there is external harmonic or 
Coulomb-like iV-body potential as given by eq. ( |33| ) and in the case of 
both An-i and BCn root systems, it is not clear if this is so in our 
case even though it is likely that this may be so in the A^-i case. 

4. Whereas our system, both on a line and on a circle, has good thermo- 
dynamic limit (i.e. E/N is finite for large N), CSM does not have good 
thermodynamic limit in either case and E/N diverges like N for large 
N. 

5. In both the cases, the norm of the ground state wavefunction can be 
mapped to the joint probability density function of the eigenvalues of 
some random matrix. Using this correspondence, in both the cases, one 
is able to calculate one- and two-point functions. However, whereas in 
the CSM this is possible only at three values of the coupling (corre- 
sponding to orthogonal, unitary or simplictic random matrices), in our 
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case the correlation functions can be computed analytically for any in- 
tegral or half-integral values of the coupling while numerically it can 
be done for any positive (3. 



6. In the CSM case with an external potential of the form 

v (e* 2 ) = ^X>? + b (e^ 2 ) 2 + c (E- 



(144) 



it has been shown [24| that the norm of the ground state wave func- 



tion can be mapped to a random matrix corresponding to branched 
polymers. It is not known if a similar mapping is possible in our case. 

7. A multi-species generalization of CSM has been done [Q, it is not clear 



if a similar generalization is possible in our case or not. 

Generalization to D-dimensions {D > 1) is possible in CSM as well 
as in our model and in both the cases one is able to obtain only a 
partial spectrum including the ground state. In both the cases, both 
two- and three-body interactions are required. Whereas our system 
has a good thermodynamic limit in any dimension D, the CSM does 
not have a good thermodynamic limit in any dimension. However, 
whereas the norm of the ground state wave function can be mapped to 
complex random matrices in the CSM case in two dimensions [|3l[], no 
such mapping has so far been possible in our case for D > 1. 

Model with novel correlations is possible in two dimensions in both 
the cases [EJ2J but unlike CSM, our system has a good thermodynamic 



limit. 

10. In the CSM, it has been possible to obtain the entire spectrum alge- 



braically by using supersymmetry and shape invariance [34]. It would 
be nice if similar thing can also be done in our model. Further, in 
the CSM, one has also written down the supersymmetric version of the 



model [35]. It would be worth enquiring if a similar thing can also be 



done in our model. 

11. In the CSM type models, one knows the various exactly solvable prob- 
lems in which the iV-particles interact pairwise by two body interaction 
36|j . The question one would like to ask in our context is: what are 



the various exactly solvable problems in one dimension in which the N 
particles have only nearest- and next-to-nearest neighbour interactions? 
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12. In the CSM, not only one- and two-point but even n-point correlation 
functions are known. It would be nice if the same is also possible in 
the present context. 



13. A la Haldane-Shastry spin models [p?[| , can we also construct spin mod- 
els in the context of our model? 

14. Unlike CSM, in our case the off-diagonal long-range order is nonzero in 
the bosonic version of the many-body theory in one dimension. Note 
however that the off-diagonal long-range order is nonzero in the CSM 
in two dimensions. 

Appendix 

1. Proof of the representation ( |84]) of Cn 

By construction, the square of the wave function (|69| ) is a symmetrical 
function of all its arguments, so that we can write eq. (|74"D as well: 

■(x u ...,x N ) | 2 , (145) 

where the particle coordinates are now properly ordered. We are thus 
allowed to substitute 4>n for ip N in ( |145| ) and obtain from eqs. (ffOl) and 







/ dx\ I dX2-- 


/ dx N 


Jo Jo 


Jo 



f 1 f* 1 f x N-l J* 

C N = N\ / dxi / dx 2 ... / dx N S{x n - x n+ i) . (146) 
Jo Jo Jo 

Changing the integration variables (xi, xn) to (Ai, A 2 , Ajv-i, xn) 
where 

An = x n -x n+ i] (n = 1, ...,N - 1) , (147) 
one easily gets 

f l rl-Ai /4-Ai-...-Ajv_a 

C N = N\ dA 1 / dA 2 ... / dAjv-l X 

Jo Jo Jo 

iJV— 1 



X / ^ P=1 dxNllSiA^SixN-x,) 2 . (148) 

J ° n=l 
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Since xn — x\ = — J2p=i A p is in fact independent of xn in the new 
set of variables, eq. ( |148|) becomes, using also S(—x) — S(l — x): 

pi rl-Ai / -l-Ai-...-A JV _ 2 

C N = N\ dA 1 / dA 2 ... / dA N ^ X 

Jo Jo Jo 

N-l N-l N-l 

X (l-EA P )Il^An) 2 5(l-EA P ) 2 . (149) 

p=l n=l p=l 

It is now convenient to introduce the extra variable 

JV-l 

A^v = 1-E A p> ( 15 °) 
P =i 

and to recast eq. ( |149| ) in the form 

pi pi pi pi N 

C N = N\ dA 1 / dA 2 ... / rfA tf _! / dA N 5(l - V A p ) X 
Jo Jo Jo Jo ~J 



N 



n=l 

AT i JV A 7 



X Aat JJ S(Aj 

n=l 

= N\ / dAi... / dA N 5{l - E A p )- £ A m J] 5(A„) i 

1/0 1/0 p=l ^ m=l n=l 

11 A 7 AT 

= (iV-1)!/ dA x ... / dA N 5(l-J2\)U S ( A n)\l51) 

J ° J ° p=l n=l 

In the second equality, we have used the fact that, apart from the factor 
An, the integrand and the integration range are completely symmet- 
rical in the variables (Aj., Ajy). Finally, the integration over these 
variables factorizes after introducing the representation 

5(1 - £ A p ) = -L r dxe -^-Y,U ^ , (152) 

p=l ^7T J-00 



and interchanging the x- and A— integrations. This produces eq. (p4|) . 

2. Proof of the representation (|§5|) of Ajy 

Proceeding along the same lines, we first put the expression flgip of An 
in the form 

A N = (N — 1)! / rfxi / dx 2 ... I dx N _i(j) N (x 1 , x N -i, 0) A 

JO JO ■/ 

X i2jv(xi, ...,xjv_i) , (153) 
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where 



R N (xi, ...,a;jv_i) = / dx(j> N (xi, ...,x N -. 1: x) 
Jo 

±/ dx(j) N (xi, ...,x, x N -i) + ... + / dx(j) N {x : xi, ...,xjv-i) 

Jxpf—1 J XI 

dx(j) N (x 1 , ...,x N - U x) + / dx<j) N (x,xi,...,x N -i) 

N - 2 rx p 

J2 u p dx(/) N (xi, x p , x, x p+1 , x N ) . (154) 

p=l "'a-V+i 



Here, i/p = l(z/ p = (— l) p ) for bosons (fermions) and we have used the 
restriction to odd iV in the second case. Thanks to the periodicity and 
the cyclic symmetry of <pN > the first two terms in the last expression 
above can be collected to give 



XJV-1 



dxcj) N (x,x 1 , ...^jv-i) 
Hence Rn becomes (with x^ = x\ — 1) 



N-l 



i^jv (zi, ...,xjv_i) = z2v p dx<p N (xi, ...,x p ,x,Xp+i, ...,xn-i) 

p=\ Jxp+i 
N-l N-l Xp 

= X) v p II ~~ / - Xp+i)^^ - x) , (n^p) 

p=l n =l Jxp+\ 
N-l N-l 

= X v v X S ( x n ~ x n+1 )S 2 (x p - x p+ i) , (n^p) (155) 
p=i n=i 

according to the definition (p3[). We also have: 



AT-2 



(xi, ...,xjv-i,0) = fl S(x m - x m+ i)S(x N _i)S(xi) . (156) 



m=l 



Inserting eqs. ( |155| ) and ( 156|) in eq. Q153 ) and introducing as before 
the new integration variables A n = x n — x n+ \ (n — 1, N — 2) and 
2at_i, we obtain 

rl /-1-Ai /-I-Ai-.-.-Ajv-s 

A* = (N-l)\ dA x dA 2 ... dA N _ 2 X 

Jo Jo Jo 
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rA N ^ N-2 

X / dx N _i [[ S , (A m )S , (x A r_i)S'(xAf_i - Aat_i) X 

J ° m=l 



N-l N-l 

X E^IIWIM; (n^p) (157) 

p=l n=l 



where Ajv_! = 1 — J2 p =i A p . The integration over x^-i gives the 
factor S 2 (A N _i) in place of S(xN-i)S(x]y^i — A^-i), so that 

1 i N-l rN-2 

A N = (N-1)\ dAt... dA N ^5(l-J2 A p) II x 
Jo Jo ~\ Ljt_i 



p=i 

N-l N-l 



m=l 



X S(A m )S 2 (A N ^) ]>> P II S(A n )S 2 (A p 

p=l n=l 



(n^pjjl58) 



On taking into account the complete symmetry of the integration mea- 
sure, one finds that the square bracket in eq. ( |158| ) can be replaced 
by 

N-2 N-3 

[...] = n S(A m ) 2 S 2 (A N ^) 2 + m J] S(A m f X 

m=l m=l 

X [S(A N ^ 2 )S 2 (A N ^)][S(A N ^)S 2 (A N - 1 )], (159) 



where r/j^ is as defined in eq. fl87|). Finally, one obtains the factorization 
of the multiple integral in eq. ( |158| ) by using again the representation 
( |1 52| ) of the S measure (with (N — 1) in place of N). This entails eq. 



A last remark may be in order. Alternative, equivalent forms of the 
representations ( |8^] ) and fl35|) would be obtained by relying on Fourier ex- 
pansions instead of Fourier integrals, that is by considering the integrands in 
eqs. ( |149| ) and ( |157| ) not as functions with compact supports [0, 1]^ C 1Z N , 
resp. [0, l]^" 1 C TZ 1 ^^ 1 , but as periodic functions (this would amount to 
modifying eq. (|152[) accordingly). It turns out however that the resulting 



representations of Cn and (as Fourier series) are much less convenient 
for the explicit or asymptotic evaluations of these quantities. 
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Figure Legends 



Fig. 1 The two-point correlation function for four integer values of f3 (from 
left to rightmost are increasing values from 1 to 4) shows clearly an absence 
of long-range order. 

Fig. 2 The two-point correlation function for some fractional values of (3 
plotted alongwith (3 equal to 1 and 2. /,From left to rightmost are increasing 
values from 1, 4/3, 3/2, 5/3, 2, 7/3, and 5/2. Thus, even for fractional values, 
there is no long-range order. 
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